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ABSTRACT 

Blue horizontal-branch (BHB) stars with T e ff approximately larger than 1 1500 K show several 
observational anomalies. In globular clusters, they exhibit low rotational velocities, abundance 
anomalies (as compared to cluster abundances), photometric jumps and gaps and spectro- 
scopic gravities lower than predicted by canonical models. It is commonly believed that the 
low rotational velocities of these stars permit atomic diffusion to be efficient in their atmo- 
sphere thereby causing the observed anomalies. Recent detections of vertical stratification of 
iron (and some other chemical elements) in several BHB stars concur with this framework. In 
this paper, improved model atmospheres that include the vertical stratification of the elements 
are applied to BHB stars to verify if they can explain their observational anomalies. The re- 
sults from theoretical model atmospheres are consistent with the photometric jumps and gaps 
observed for BHB stars in globular clusters. It is found that iron stratification in the theo- 
retical models and that obtained from observations have similar tendancies. Our results also 
show that the spectroscopic gravities obtained while using chemically homogeneous model 
atmospheres to fit observations are underestimated. These results significantly strengthen the 
belief that atomic diffusion is responsible for these BHB-star anomalies. 
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1 INTRODUCTION 

Horizontal-branch stars are evolved intermediate-mass stars that 

r* are burning helium in their core (e.g. Hoyle & Schwarzschild 

. J^ 1955). Some hot blue horizontal-branch (hereafter BHB) stars with 

K^ r e ff larger than approximately 1 1500 K are very interesting objects 

$_^ since they exhibit several observational anomalies. This paper aims 

Cw to study certain aspects of these astronomical objects in the light of 

recent modelling and observational results. 

First, abundance anomalies are observed for such BHB stars 
in several globular clusters (Glaspey et al. 1989; Behr et al. 1999; 
Moehler et al. 1999; Behr et al. 2000a; Behr 2003a; Hubrig et al. 
2009). Khalack et al. (2007, 2008 and 2010) have detected vertical 
stratification of the abundance of several chemical elements includ- 
ing Fe in the atmosphere of some stars of this type. Khalack et al. 
(2008) found that the iron abundance increases toward the lower 
atmosphere in three BHB stars. Nitrogen and sulfur stratification 
was also detected in the hot BHB star HD 135485 (Khalack et al. 
2007). The abundances of these two elements are found to increase 
toward the upper atmosphere. More recently, Khalack et al. (2010) 
determined the vertical Fe abundance gradient in a sample of 14 
BHB stars using the observed spectra of Behr (2003a). They de- 
tected vertical Fe stratification in five (and possibly seven) of these 
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BHB stars. Isotopic anomalies have also been recently detected in 
BHB stars (Hubrig et al. 2009). 

Photometric jumps are observed on the hot side of 7" e ff ~ 
11500 K for the horizontal-branch sequence of several globular 
clusters (Grundahl et al. 1999) as compared to what is predicted 
by canonical models. Photometric gaps are also observed for this 
sequence (Ferraro et al. 1998) at the same T e g where the photomet- 
ric jump mentioned above occurs. 

Another anomaly detected is that the rotational velocities of 
BHB stars are observed to drop abruptly for T^ approximately 
larger than 11500 K (Peterson, Rood & Crocker 1995; Behr et al. 
2000a and 2000b; Behr 2003b). 

Finally, spectroscopic gravities of BHB stars in globular clus- 
ters, at least for the metal-poor ones (e.g. Crocker et al. 1988; 
Moehler et al. 1995), are also lower than those predicted by classi- 
cal BHB models. 

Atomic diffusion (Michaud 1970) is effective only if the 
medium is hydrodynamically stable enough to prevent mixing, be- 
cause of the order of magnitude of the diffusion velocities which 
are much smaller than those of macroscopic motions. In the case 
of BHB stars, their relatively low rotational velocities should ren- 
der their superficial layer stable enough for diffusion to take place. 
For instance, recently, Quievy et al. (2009) showed that for val- 
ues of rotational velocities found for BHB stars with T e ff above 
1 1500 K, meridional circulation in not efficient enough to prevent 
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He from gravitationally settling. This leads to the disappearance 
of the superficial He convection zone and renders the atmosphere 
more stable. Once the medium is stable enough, atomic diffusion 
can dominate there, leading to abundance anomalies and vertical 
abundance stratification of the chemical elements through the mi- 
gration of the particles caused by the disbalance of gravity and ra- 
diative forces, and yielding superficial abundance anomalies. Strat- 
ification of the elements has an effect on the physical structure of 
the atmosphere through the change in the opacities, thus leading to 
photometric anomalies, and to a change of the shape of the Balmer 
lines, explaining the low values of surface gravity obtained by fit- 
ting of these lines, if one uses canonical models. The scenario out- 
lined above, where atomic diffusion takes place in the atmospheres 
of hot BHB stars, may therefore explain the various observational 
anomalies observed for these stars. 

Hui-Bon-Hoa, LeBlanc and Hauschildt (2000) constructed 
stellar atmosphere models of BHB stars with vertical stratification 
of the elements, and self-consistent atmospheric structure. These 
models were successful in qualitatively reproducing the above- 
mentioned anomalies when assuming that atomic diffusion be- 
comes efficient in BHB stars with T e ff > 11500 K (i.e. the tem- 
perature above which the vast majority of the BHB stars rotate 
slowly). These models have been improved recently (LeBlanc et 
al. 2009), and the aim of this paper is to apply these new mod- 
els to BHB stars and to compare the results to observational data. 
Among the improvements, the convergence scheme to obtain self- 
consistent abundance stratification was ameliorated in these mod- 
els. Also, the diffusion coefficients now used take into account the 
interaction between the ions and neutral hydrogen. Moreover, in 
these new models, upper and lower abundance limits were imposed 
to avoid that certain elements attain extremely large (for the test- 
particle approximation of the diffusion theory to still be valid) or 
small (to avoid numerical problems) abundances at certain depths 
(see Section 3 of LeBlanc et al. 2009 for more details). 

First the new model atmospheres with elemental stratification 
will be briefly described. Theoretical vertical stratification of Fe in 
models of various T e ff will be compared to the observational re- 
sults of Khalack et al. (2007, 2008 and 2010). The synthetic pho- 
tometry obtained with these models will then be presented, check- 
ing if they are able to reproduce observed photometric jumps and 
gaps for globular cluster BHB sequences. Synthetic Balmer lines 
are also compared between models with different prescriptions, to 
show how the diffusion model can account for the low values of 
gravity derived when using canonical models. 



2 MODEL ATMOSPHERES WITH STRATIFIED 
ABUNDANCES 

The process of elemental stratification relies on the diffusion ve- 
locity, which causes the migration of the different chemical ele- 
ments within stars. The leading terms that cause this velocity are 
gravity and the radiative acceleration resulting from the momentum 
transfer between the radiation field and each chemical species (e.g. 
Gonzalez et al. 1995). This transfer depends on the opacity of the 
species under consideration and the local monochromatic radiation 
field, which, in turn, via the monochromatic opacities, depends on 
the local abundances of the different species. The stratification of 
the abundances is thus a time-dependent process and it is therefore 
complex to treat it thoroughly. As the radiative transfer equation 
has to be solved explicitly in stellar atmospheres, which are op- 
tically thin media, time-dependent calculations are not feasible at 



Table 1. Physical parameters of models. 
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M(M, S ) 


11000 


3.78 
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12000 


3.95 
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4.11 


0.625 
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4.26 


0.609 


16000 


4.52 


0.581 


18000 


4.74 


0.561 



the present time. Studies of abundance stratification can however 
be performed assuming that an equilibrium state can be reached, 
when the diffusion velocity of each element is nil at each layer 
of the atmosphere. This is the framework within which our cal- 
culations are performed. Even though the elemental stratification 
profiles can possibly differ from those computed time-dependently, 
this approximation can be used to gauge the impact of vertical strat- 
ification of the abundances on the atmospheric structure and on cer- 
tain observed quantities. For a more detailed discussion surround- 
ing time-dependent diffusion and equilibrium solution, the reader 
is referred to Alecian & Stift (2007) in which Ap stars are studied. 

The model atomspheres presented here are calculated with 
a modified version of the PHOENIX code (Hauschildt, Allard & 
Baron 1999) as described by LeBlanc et al. (2009). The models 
computed here are in LTE and include 39 elements (H-Ga, Kr-Nb, 
Ba and La). To build the stratification profile for each element, the 
code seeks iteratively a solution where the atmospheric structure 
and the vertical abundance stratifications yield a nil diffusion ve- 
locity for each chemical species (i.e. the equilibrium solution). In- 
deed, as the abundance of each species changes locally, so do the 
opacities, and thus the atmospheric structure. The resulting model 
atmospheres are thus self-consistent from this point of view. 

For the evaluation of the radiative accelerations, both bound- 
bound and bound-free transitions are considered for each species 
while assuming LTE (Hui-Bon-Hoa et al. 2002). As the radiative 
transfer is solved with the opacity sampling method, the radiative 
accelerations, along with the radiative flux, are calculated using a 
sufficiently fine frequency grid to ensure adequate precision (e.g. 
LeBlanc, Michaud & Richer 2000). Several improvements were 
brought to the models presented here as compared to the models 
used in Hui-Bon-Hoa, LeBlanc & Hauschildt (2000). For example, 
the convergence scheme used to compute the stratification of the 
abundances has been modified. Also, more precise diffusion coeffi- 
cients are used. More details concerning the improvements brought 
to the atmospheric code are found in LeBlanc et al. (2009). 

A series of BHB model atmospheres from T e ff = 11000 to 
18000 K (see Tab.fTb were constructed while including abundance 
stratifications. The surface gravities and masses used here are those 
of the zero-age horizontal-branch (ZAHB) models of VandenBerg 
et al. (2000). 

The stratification of the elements predicted at equilibrium 
strongly modifies the physical structure of the stellar atmospheres 
of hot BHB stars. Figure [T] shows the ratio of the temperature as 
a function of optical depth in several of these models as compared 
to the temperature of corresponding models (with the same funda- 
mental parameters) with homogeneous models with a metallicity of 
-1.5 dex solar. The abundances used for these homogeneous models 
are chosen to be -1.5 dex compared to solar abundances which are 
typical of well-studied globular cluster metallicity (e.g. M13, M3 
etc.). The large modifications in the structure of the self-consistent 
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Figure 1. The ratio of the temperature as a function of optical depth in 
various models including stratification to that of a chemically homogeneous 
model with metallicity of -1.5 dex solar of the same effective temperature 
and surface gravity. The dotted-dashed line shows the position of a ratio of 
one. 



Figure 2. Abundance of Fe relative to optical depth for self-consistent 
model atmospheres of BHB stars with T af from 1 1000 K to 18000 K. 



model atmospheres in the line-formation region as compared to ho- 
mogeneous ones can give rise to detectable observational effects. 
As will be discussed below, the diffusion process not only causes 
abundance stratification of the elements but can also lead to visible 
effects related to the photometric colors and spectroscopic gravi- 
ties of hot BHB stars. The results of these theoretical models are 
compared below to well-established anomalies of hot BHB stars. 



3 STRATIFICATION OF IRON IN BHB STARS 

The vertical stratification of Fe in our stratified model atmospheres 
is shown in Fig. [2] The results presented in this figure clearly show 
that the gradient of the Fe abundance (as a function of optical depth 
in the atmosphere) in the range -4 < log T5000 ^ "2 decreases with 
the increase of T e s. Since this range of optical depth is where a lot 
of the iron lines are formed, such a tendancy can be verified with 
spectroscopic studies. 

Figure 15] compares the slopes of the detected Fe abundance 
gradients or lack thereof for a large number of BHB stars (Khalack 
et al. 2007, 2008 and 2010) to those predicted by the theoretical 
model atmospheres of LeBlanc et al. (2009). In these observational 
studies, the vertical stratification of iron was gauged by determining 
the abundance of individual lines formed at various optical depths. 
The dots in this figure represent individual hot BHB stars while the 
solid line is a quadratic fit of the observed abundance slopes. The 
filled diamonds connected by a dashed line give the linear slopes 
of Fe abundance in the theoretical models. These slopes are those 
of the linear fit of the abundance stratification of iron for the at- 
mospheric layers in the range -5 < log X5000 < -2 for each T e ff. 
These layers correspond to the atmospheric depths where most of 
the Fe lines are formed and where iron stratification is detected (see 
Figure 1 of Khalack et al. 2010). Figure |3j shows that the observed 
slopes of the Fe abundance decreases with T e ff for hot BHB stars 
and eventually becomes undectectable at T e g ~ 14000 K, in agree- 
ment with the theoretical results of the models of LeBlanc et al. 
(2009) shown in Fig. [2] These results may serve as additional proof 
that elemental stratification due to atomic diffusion occurs in the 
atmospheres of hot BHB stars. 
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Figure 3. Slope of the vertical Fe abundance obtained from the analysis 
of observed spectra (filled dots with their corresponding error bars) taken 
from Khalack et al. (2007, 2008 and 2010). Also shown in this figure is a 
quadratic fit of these observed slopes (solid line) and the linear slopes (filled 
diamonds joined by dashed line) obtained in the theoretical models in the 
range -5 < log T5000 < -2 which corresponds to the depth from where most 
of the Fe lines emanate. The dotted line shows the position of a nil slope. 



4 PHOTOMETRIC JUMPS AND GAPS 

BHB stars in globular clusters exhibit several photometric anoma- 
lies namely photometric jumps and gaps. This section aims to in- 
terpret them in the light of self-consistent stratified model atmo- 
spheres. 

Grundahl et al. (1999) observed a photometric jump in the 
(u,u — y) color-magnitude diagram for the horizontal branch se- 
quence in several globular clusters. This jump consists of a shift 
towards brighter stars of the horizontal-branch sequence with re- 
spect to canonical models, and occurs for stars with T e ff higher than 
approximately 11500 K. Therefore, there exists a difference (for 
r e ff y 11500 K) between the observed and predicted colors. The 
PHOENIX model atmospheres including elemental stratification of 
Hui-Bon-Hoa, LeBlanc & Hauschildt (2000) were able to qualita- 
tively explain this jump, when stars with T e g larger than 11500 K 
become stable enough to allow atomic diffusion to have visible ef- 
fects. This section aims to verify if the new models of LeBlanc et 
al. (2009) also reproduce this jump. 
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Figure 4. Color-magnitude diagram («, u — y) for the several sets of BHB- 
star models. The filled dots represent models from T e ff = 1 1000 to 18000 K 
described in Tab. [T](r e ff increases from right to left in this figure) that in- 
clude vertical elemental stratification while the empty dots represent mod- 
els with homogeneous abundances (with metals at - 1 .5 dex their solar val- 
ues). An arrow shows the photometric jump predicted if diffusion becomes 
efficient at T c ff = 11500 K. The dotted curve represents atmosphere with 
homogeneous abundances equal to solar values. The dashed curve are the 
results when assuming stratification but without changing the structure of 
the atmosphere. 
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Figure 5. Color-magnitude diagram (V, V — V) that shows the presence of a 
gap assuming that diffusion becomes efficient at 7" e ff = 1 1 500 K. The curves 
are defined in the caption of Fig. [4] 



FigureEJshows the («, U—y) color-magnitude diagram for stars 
with and without elemental stratification. These results are consis- 
tent with the observational results of Grundahl et al. (1999) when 
compared to canonical models. Our results show that if diffusion 
becomes efficient at T e ff ~ 1 1 500 K, a jump occurs in the (u,u — y) 
diagram. The difference in the u-bandpass magnitude between the 
homogeneous models (at -1.5 dex solar abundance for the metals) 
and those with stratification is 0.5 mag or more. 

Two other curves are also shown in this figure. The results 
while assuming that there is general enrichment of the metals (as 
compared to average cluster abundances) equal to a solar value is 
shown (dotted curve). This curve clearly shows that a general en- 
richment of the metals throughout the atmosphere cannot account 



for the photometric jump observed. The other curve shown in Fig- 
ure H] (dashed curve) represents the colors when one calculates the 
stratification profiles in the stars but without modifying the atmo- 
spheric structure (i.e. keeping the structure equal to that of the ho- 
mogeneous model with -1.5 dex solar abundances). This shows 
that opacity changes cannot by themselves explain the photomet- 
ric jump for BHB stars and that the structural changes are essential 
for this anomaly to occur. In summary, the photometric changes of 
the stratified models as compared to the homogeneous models are 
brought about by the change of the atmospheric structure due to 
the vertical stratification of the elements. For instance, the temper- 
ature in the line forming regions of hot BHB stars (see Figure fT} 
found in the models with stratification is significantly larger (by up 
to approximately 20%) as compared to the models with an homo- 
geneous composition with -1.5 dex metal solar abundances. 

Photometric gaps are also observed on the horizontal-branch 
sequence of globular clusters. These gaps in color-magnitude di- 
agrams are characterized by an underpopulation of stars for spe- 
cific values of colors. For example, a photometric gap is observed 
in the (V, U — V) color-magnitude diagram (Ferraro et al. 1998). 
This gap corresponds to T & g ~ 1 1500 K where rotation plumets for 
BHB stars (Peterson, Rood & Crocker 1995; Behr et al. 2000a and 
2000b; Behr 2003b) and where diffusion could begin to be efficient. 
Hui-Bon-Hoa, LeBlanc and Hauschildt (2000) had shown that this 
gap may be explained if diffusion becomes efficient at T e ff ~ 1 1500 
K. The new models presented here also show the same feature (see 
Fig.[5](. The width of the gap is more than 0. 1 mag which is consis- 
tent with the observations (see gap named Gl in Fig. 3 of Ferraro 
et al. 1998). As in FigureH] two other curves are plotted as well to 
show that the photometric gap cannot be explained by a homoge- 
neous enrichment of the metals nor elemental stratification without 
structural changes. 



5 SPECTROSCOPIC GRAVITIES 

Several spectroscopic studies of BHB stars in metal-poor globular 
clusters give surface gravities lower than those predicted by canon- 
ical ZAHB models for stars with temperatures between 1 1,000 and 
20,000 K (see Crocker et al. 1988; Moehler et al. 1995). The values 
of the surface gravities are determined through the study of Balmer 
lines profiles. 

To see the effect of stratified model atmospheres on the value 
of surface gravities obtained this way, we synthesized Balmer lines 
for the self-consistent model atmospheres with the different effec- 
tive temperatures used in the previous sections. Figure [6] summa- 
rizes our results. Computed with values of surface gravity of canon- 
ical horizontal-branch models (see Tab. [lb, the stratified models 
show Balmer line profiles closer to those of homogenous models 
of the same effective temperature but with lower gravity (by ap- 
proximately 0.4 to 0.5 dex, see dashed line of Figure[6l. For exam- 
ple, for the 14000 K model, the homogeneous model with log g = 
3.8 (dashed line) better fits the wings of the Hr line predicted by 
the self-consistent stratified model than the homogeneous model 
with the canonical value of log g = 4.26 (dotted line). Thus, try- 
ing to determine spectroscopic gravities using homogenous mod- 
els is misleading when the star has its structure altered by abun- 
dance stratification. This effect mimics lower spectroscopic gravi- 
ties when using homogeneous models as a workbench. The use of 
model atmospheres taking into account abundance stratifications is 
therefore critical for spectroscopic gravity determinations. 

As a test (not shown in Figure 16), we also computed Balmer 
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Figure 6. Hp line profiles for the 1 1000 and 14000 K models. The stratified 
model profiles are denoted "stratified" and assume repectively log g = 3.78 
and 4.26 (see Tab.fTJ, whereas homogenous cluster abundance models (-1.5 
dex solar) are named following their value of log g used. The first value 
(dotted line) is that of the canonical sequence, and the second (dashed line) 
that of the homogenous model that fits best the stratified model. 



ical models predict that the iron abundance stratification decreases 
as a function of r en - and becomes negligible for BHB stars hotter 
than T e ff ~ 14000 K. This result is consistent with the detected iron 
abundance slopes found by spectral line analysis of BHB stars un- 
dertaken by Khalack et al. (2007, 2008 and 2010). 

The model atmospheres with vertical stratifications also repro- 
duce the photometric jumps and gaps that occur for BHB stars at 
r e ff ~ 1 1500 K. These anomalies were the precursors for the belief 
that abundance stratification due to diffusion is present in the slowly 
rotating hot BHB stars. The theoreticals models presented here also 
predict, as previously observed, that the spectroscopic gravities ob- 
tained for these stars while using homogeneous model atmospheres 
to fit observations systematically underestimate their true surface 
gravity. The photometric and surface gravity anomalies are due to 
the large modification in the atmospheric structure brought about 
by the stratification of the elements in the atmosphere. 

In summary, the results presented here reaffirm those found by 
Hui-Bon-Hoa, LeBlanc & Hauschildt (2000) regarding the photo- 
metric jumps and gaps as well as for the lower spectroscopic grav- 
ities found in hot BHB stars. In addition, we have found that the 
general trend for the observed Fe stratification slopes in hot BHB 
stars is similar to the one predicted by the theoretical models (i.e. 
the slope of Fe abundance decreases as T e ff increases). 

The theoretical models used here assume that the elements 
attain an equilibrium where their respective diffusion velocities 
tend towards zero. The approach used to calculate the model at- 
mospheres is not a time-dependent treatment of the diffusion phe- 
nomena. Such time-dependent modelling of the atmosphere is es- 
sential to obtain more precise elemental stratification profiles. Since 
Michaud et al. (2008) have calculated evolutionary models for BHB 
stars while including atomic diffusion, their results could be useful 
to constrain the amount of each element exchanged with the stellar 
interior, when building up time-dependent elemental stratifications 
in the atmospheres. 

Even with the uncertainties inherent in the treatment of the 
diffusion phenomenon in the theoretical models used here, the re- 
sults presented in this paper significantly strengthens the belief that 
atomic diffusion is responsible for the aformentioned anomalies 
observed for hot BHB stars. In the near future, we plan to inves- 
tigate the error engendered by the use of homogeneous models for 
the photometric determination of fundamental parameters of BHB 
stars. Such a study is critical for ensuring the use of proper funda- 
mental parameters for this type of star. 
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